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Section  1 


INTRODUCTION 

Under  the  sponsorship  of  the  Meteorology  Oivision  of  the  Air  Force 
Geophysics  Laboratory  (AFGL)  we  have  undertaken  the  development  of  a 
parameterization  program  involving  the  transfer  of  thermal  infrared  and 
solar  radiation  in  clear  and  cloudy  atmospheres.  The  prime  objective 
of  this  program  is  to  construct  a  unified  and  coherent  radiation  package, 
including  cloud  contributions,  which  is  computationally  economic  yet 
numerically  reliable  and  stable,  for  the  incorporation  into  a  global 
dynamic  model.  Thus,  our  radiation  and  cloud  programs  are  closely  in 
line  with  the  overall  AFGL  effort  in  developing  a  global  model  for  short 
and  medium  range  weather  and  cloud  forecast  research. 

In  this  report,  we  describe  the  infrared  and  solar  radiation  param¬ 
eterization  programs  developed  at  the  University  of  Utah  over  the  period 
supported  by  AFGL.  In  Sections  2  and  3,  we  present  detailed  descriptions 
of  the  physical  concepts  and  basic  parameterization  equations  developed 
for  the  computation  of  infrared  and  solar  fluxes  in  clear  and  cloudy 
atmospheres.  Verifications  of  the  IR  and  solar  radiation  parameteriza¬ 
tion  programs  also  are  given  In  these  sections.  In  addition,  we  present 
in  Section  4  some  preliminary  results  of  the  global  distributions  of 
clouds  and  radiation  budgets  generated  from  the  AFGL  model  utilizing 
the  atmospheric  data  corresponding  to  15  January  1978.  Finally,  detailed 
documentation  of  the  radiation  parameterization  programs  including  the 
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programming  Interface  with  the  AF6L  six-layer  general  circulation  model 
is  given  in  the  Appendix. 


Section  2 

INFRARED  RADIATION  PROGRAM 
2.1  Clear  Atmosphere 

Following  the  analysis  presented  by  Liou  (1980)  utilizing  the  broad¬ 
band  emissivity  concept,  the  upward  and  downward  infrared  fluxes  at  a 
given  level  in  a  clear  atmosphere  may  be  written  in  the  forms 

Ff(z)  =  oT4  [1  -  ef(z,T*)]  +  f  oT4(z' )  K(z-z')  dz1  (2.1) 

F*(z)  =  f2  oT4(z' )  K(z'-z)  dz*  (2.2) 

Jzco 

where  the  broadband  flux  emissivity  and  the  atmospheric  kernel  function 
are  defined  by 

f  6  f 

e  (z,T*)  -  T  ej  [a1(z),T,)3  (2.3) 

K(z-z')  »  -  j|r  €f  [z-z',T(z’)] 

=  -  gfz  f  ej  [u,(z)  -  a,(z'),T(z')]  (2.4) 

In  these  equations,  T*  denotes  the  surface  temperature,  uj  “  u2  *  u3  = 
u  ,  the  water  vapor  path  length,  u*  3  u  ,  the  carbon  dioxide  path  length, 

W  *T  W 

and  Uj  3  uQ,  the  ozone  path  length,  (i  3  1-5)  are  the  individual 
broadband  emissivities  for  the  HgO  rotational  band,  HgO  vibrational - 
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rotational  band,  H20  continuum,  C02  and  03  bands,  respectively.  The 
broadband  emissivity  corresponding  to  the  H20-C02  overlap  correction 
is  given  by  e^Ug)  =  -e[(uw)  ej(uc).  In  Eqs.  (2.3)  and  (2.4),  Q.  denotes 
the  pressure  and  temperature  corrected  path  length  defined  below  for 
each  band.  The  individual  broadband  emissivities  are  defined  by 


4(0,, T)  -  l  hB-  (T)  [1  -  4(a  )]Av/o  T4  , 
n  n  «  1 


(2.5) 


where  n  denotes  the  number  of  spectral  intervals  Av  ,  irB-  the  spectral 

n  r 

Planck  flux,  and  the  spectral  slab  (or  diffuse)  transmittance  T^(u)  a  t- 
(1.66  u).  The  spectral  transmittance  T-  is  obtained  either  from 
theoretical  band  models  or  from  detailed  absorption  line  data  involving 
the  line  position,  line  strength  and  half  width. 


2.1.1  Water 


For  the  three  water  vapor  absorption  bands  we  use  the  broadband 
emissivity  curves  in  terms  of  polynomial  functions  presented  by  Liou 
and  Ou  (1983).  In  their  approach,  the  broadband  emissivities  for  water 
vapor  pure  rotational  and  vibrational -rotational  bands  are  derived 
from  the  random  band  parameters  assuming  Lorentz  line  shapes  obtained 
by  Rodgers  and  Walshaw  (1966)  in  which  an  empirical  temperature  cor¬ 
rection  for  the  rotational  band  was  also  given.  For  the  water  vapor 
continuum  absorption  in  the  window  regions,  the  empirical  formula 
developed  by  Roberts  et  al .  (1976)  was  used.  The  polynomial  functions 
have  the  form 


In  el(Q.,T)  =  )  a„  .  u1?,  u.  =  (2  log,„  0.  -  a)  /  b 


(2.6) 


STOVE 


10  -2  -7  -2 

with  u.  max  (=  10  g  cm  )  and  u.  m^n  (=  10  g  cm  )  representing  the 

maximum  and  minimum  path  lengths  used  in  the  numerical  fitting.  The 

pressure  and  temperature  corrected  path  lengths  are  (Liou  and  Ou,  1981) 


where  A'  =  10.6278  x  103  K1 ,  B'  «  -44.6132  x  10'6  K'2,  T  =  260°  K, 

a 

Tb  =  296°  K,  Pq  =  1013  mb,  Pw  and  P^  are  vapor  pressure  at  temperatures 
T  and  Tfa,  respectively  and  pw  is  the  water  vapor  density.  Table  1  lists 
the  coefficients  of  the  broadband  emissivity  equation  an  .  for  three  water 
vapor  bands  (i  =  1,2,3).  As  illustrated  in  the  paper  by  Liou  and  Ou  (1983), 
the  clear  column  cooling  rates  due  to  water  vapor  absorption  and  emis¬ 
sion  in  the  troposphere  agree  well  with  results  from  line-by-line  calculations 
Pertinent  results  will  be  duplicated  in  this  report  for  documentation  purposes 

2.1.2  Carbon  Dioxide 

The  broadband  emissivity  for  carbon  dioxide  is  computed  using  the 
parameterization  scheme  recently  developed  by  Ou  and  Liou  (1983).  In 
that  paper  the  COg  broadband  emissivity,  which  is  a  function  of  the 
pressure  and  temperature  corrected  path  length  u  and  the  atmospheric 

v 

temperature,  is  derived  from  the  line-by-line  transmittance  data  pre¬ 
sented  by  Fels  and  Schwarzkopf  (1981).  It  is  fitted  by  the  following 
equations: 


& 


(2.10b) 


£n 


efo.T) 


bn  +  b. 


O’, 


for  u  < 


10"4  g 


cm 


where 


O'  *  (2  log  u  +  7.69897)  /  6.30103 
and  aQ  = -4.00893  +  f(T),  ^  =  4.39828,  a£  =  -3.07709,  a3  =  0.94529, 
bQ  =  -4.00036  +  f (T) ,  and  bj  =  5.13453.  f(T)  =  in  (h(T)  g(T)/[h(T0) 

g(T0)]>,  where  h(T)  =  1  +  A  AT(1  -  BAT),  g(T)  =  f  ‘  nBv(T)  dv/oT4,  A  =  1.833 

-4  -2  ^Vli  V  -1 

x  10  ,  B  =  1.364  x  10  ,  AT  =  T  -  250,  =  500  cm-1  and  =  850  cm  The 

pressure  and  temperature  corrected  path  lengths  have  the  form 

u4  =  2  c  /hi  +  4(c2  /  u2  +  c  /  U)]1/2  -  lj  (2.11) 

-4  -2 

where  the  constant  c  =  3.7551  x  10  g  cm  and 


Pc(z’)  dz' 


(2.12) 


with  pc  the  carbon  dioxide  density  and  Tg  =  273°  K.  CO2  cooling  rates 
calculated  from  the  broadband  emissivity  parameterization  approach  show 
a  maximum  difference  of  0.5°  C/day  when  they  are  compared  with  those 
computed  from  exact  line-by-line  integrations  utilizing  standard  and 
polar  atmospheric  temperature  profiles. 


2.1.3  Ozone 

As  for  ozone,  random  model  parameters  for  the  9.6  ym  band  derived 
by  Goldman  and  Kyle  (1968)  are  utilized  to  derive  the  03  broadband  emis¬ 
sivity  which  is  also  fitted  in  terms  of  the  polynomial  function  [Eq. 


(2.6)].  However,  the  maximum  and  minimum  path  lengths  used  in  the 
-3  -7  -2 

fitting  are  10  and  10  g  cm  ,  respectively.  Because  of  the  unce 
tainty  of  applying  pressure  and  temperature  corrections  to  the  path 
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length  in  the  cooling  rate  calculation  due  to  ozone,  no  adjustment  has 
been  made.  Table  1  also  lists  the  coefficients  of  the  broadband  emis- 
sivity  equation  for  the  ozone  band.  The  calculation  of  ozone  cooling 
rates  is  certainly  an  area  requiring  further  improvement  especially 
when  the  stratosphere  is  to  be  incorporated  in  the  dynamic  model. 


2.2  Cloudy  Atmosphere 

We  shall  consider  a  two-layered  cloud  system  having  a  semi¬ 
transparent  and  non-black  high  cloud  aloft  coupled  with  a  black  cloud 
below.  Let  the  infrared  broadband  emissivity,  transmissivity  and 
reflectivity  for  the  high  cloud  be  denoted  as  ec,  tc  and  rc,  respec¬ 
tively.  For  black  clouds,  we  let  ec  =  1  and  tC  *  rc  -  0  for  simplicity 
of  treating  the  infrared  flux  exchange.  In  a  manner  similar  to,  but 
not  exactly  the  same  as,  the  parameterization  equations  given  by  Liou 
and  Ou  (1981),  we  may  express  the  upward  and  downward  fluxes  at  a  given 
level  (Fig.  1)  above  the  high  cloud,  between  high  and  progressively 
lower  clouds,  and  below  the  lower  cloud. 

For  heights  above  the  high  cloud  top,  we  have 

F*(z)  =  [Z  <jT4(z'  )  K(z'-z)  dz'  ,  (2.13) 

OO 

Ff(z)  =  rc  jZtl  oT4(z' )  K(z'+z-2ztl)  dz' 

ZQO 

+  ^  oTt2  [1  -  ef  (zbrZt2+Z‘Ztl’Tt2)]  +  |  bl  aT4(z,)  K(z’ztl+zbl 

Zt2 

+  oT4,  [1  -  ef(z-z4.1  ,T.. )]  +  f  oT4(z’)  K(z-z')  dz'  , 


configuration  of  solar  and  1R  flux  transfer  in  an  atmosphere  containing  two  cloud  layers 


Below  the  high  cloud  base  and  above  the  low  cloud  top,  we  obtain 


Ff(z)  =  oTj2  [1  -  ef  (z-zt2,Tt2)]  +  J*  oT4(z')  K(z-z’)  dz*  ,  (2.15) 

zt2 

F+(z)  =  rcJ  aT42  [1  -  cf  (2zbi'zt2"z,Tt2^  +  °t4(z')  K(2zbl-z-z')  dz 

[_  zt2 

+  tc  j  tl  oT4(z')  Mz'-z^hi-z)  dz' 

2co 

+  <  ebl  aTbl  ^  ■  ef^zbl'z*Tbl^  +  I  °t4(2')  k(z'*z)  dz'  7  » 

Zt2  (2.16) 


ztZ<z<  zbl 


Below  the  low  cloud  base,  the  fluxes  are  related  to  two  black  surfaces 
and  are  simply  given  by 

Ff(z)  =  oT4  [1  -  cf(z,T*)]  +  [Z  oT4(z‘ )  K(z-z')  dz'  (2.17) 

Jo 

F*(z)  =  oT42  [1  -  ef(zb2-z,Tb2)]  +  j  oT4(z')  K(z-z')  dz1  (2.18) 

jzb2 

0  <  z  <  z.  o 


When  there  is  only  one  cloud  generated  in  the  model,  we  set  zb2  =  z^2  = 
zbl  and  Tb2  =  Tt2  =  Tbj.  Thus,  only  Eqs.  (2.13)  -  (2.16)  are  required 
in  the  computation.  Note  that  except  for  single  high  clouds,  all  other 
cloud  types  are  considered  to  be  black  reflectively. 

The  broadband  infrared  emissivity,  reflectivity  and  transmissivity 
for  cirrus  are  functions  of  the  vertical  ice  content  following  the 


parameterization  equation  developed  by  Liou  and  Wlttman  (1979)  in  the 
form 


5  1 

R(W)  =  l  c,  W1  (2.19) 

i=0  1 

where  R(Wj  denotes  the  flux  reflectivity,  transmissivity  or  emissivity 

of  the  cirrus  cloud,  c^  the  predictor  coefficients  listed  in  Table  2,  ' 

2  -2 

and  W  the  vertical  ice  content  In  units  of  10  g  in  . 

For  partly  cloudy  conditions,  the  infrared  cooling  rate  at  each 
model  layer  is  obtained  by  weighting  linearly  the  percentage  of  the 
fractional  cloud  cover  c  and  the  percentage  of  the  clear  portion  in 
the  form 


(2.20) 


Table  2.  Approximating  predictor  coefficients  c..  for  a  cirrus  cloud 
in  the  infrared  region. 


i 

Reflectivity 

Transmissivity 

Emissivity 

0 

0.30619E  -  01 

0.73597E  +  00 

0.28568E  +  00 

1 

0.8U34E  +  00 

-0.38162E  +  01 

0.35222E  +  01 

2 

-0.18995E  +  01 

0.83288E  +  01 

-0.76275E  +  01 

3 

0.24900E  +  01 

-0.91066E  +  01 

0.83096E  +  01 

4 

-0.15805E  +  01 

0.48765E  +  01 

-0.44405E  +  01 

5 

0.37581E  +  00 

-0.10141E  +  01 

0.92209E  +  00 

This  equation  is  also  applicable  to  solar  heating  rate  calculations  in 
partly  cloudy  conditions. 


The  radiative  heating  rate  at  level  z  is  related  to  the  divergence 
of  net  fluxes  and  is  given  by 

1  I  -  *  Wfi  4  [Ft(2)  -  F*(2)]  •  (2-21> 

where  p  is  the  air  density,  Cp  the  specific  heat  at  constant  pressure 
and  Az  the  model  layer  thickness.  The  negative  and  positive  signs  are 

applicable  to  infared  and  solar  radiation,  respectively. 

I 

| 

2.3  Verification  of  1R  Parameterization  Program 

Utilizing  the  broadband  emissivity  approach  denoted  in  Sections 
2.1  and  2.2,  we  have  carried  out  a  number  of  numerical  comparisons  and 
verifications  with  more  detailed  band-by-band  and/or  line-by-line  cal¬ 
culations  (Liou  and  Ou,  1981;  Ou  and  Liou,  1983).  Below,  we  present 
some  highlights  of  these  comparisons. 

In  Fig.  2,  we  first  illustrate  comparisons  of  infrared  cooling 
results  using  band-by-band  and  broadband  approaches  for  water  vapor 
and  ozone  components.  For  the  water  vapor  component,  the  broadband 
approach  reproduces  fairly  accurately  the  cooling  rate  below  about  20 
km  when  it  is  compared  with  the  band-by-band  approach  with  differences 
within  about  0.3°C  day"*.  The  present  results  show  reasonable  agree¬ 
ment  with  those  obtained  from  a  more  exact  line-by-line  calculation. 
Note  that  cooling  results  from  the  line-by-line  method  are  available 
only  for  heights  between  3  and  25  km.  At  above  35  km  the  broadband 
value  differs  from  the  band-by-band  result  by  as  much  as  1°C  day’1. 
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This  larger  deviation  is  in  part  due  to  the  cut-off  of  water  vapor 
concentration  at  a  height  of  45  km  where  cooling  drops  to  zero.  We 
note,  however,  that  cooling  due  to  water  vapor  at  about  20  km  is 
relatively  insignificant  compared  with  that  of  carbon  dioxide  and 
ozone.  As  for  the  ozone  component,  agreement  between  results  from 
the  broadband  and  band-by-band  methods  are  excellent  with  differences 
held  to  within  0.2°C  day"1.  Note  that  ozone  shows  heating  on  the  order 
of  1°C  day-1  at  ~25  km,  corresponding  to  the  ozone  maximum. 

In  our  previous  program  (Liou  and  Ou,  1981),  we  used  the  random 
band  parameters  calculated  by  Rodgers  and  Walshaw  (1966)  for  H20  and 
C02  to  derive  the  broadband  emissivities  and  overlap  contributions. 
Since  they  only  gave  one  C02  band  value,  the  C02  broadband  emissivity 
computed  from  one  band  value  is  generally  less  accurate,  especially 
for  atmospheres  above  about  35  km.  For  this  reason,  we  have  improved 
the  parameterization  of  the  C02  broadband  emissivity  based  on  line-by¬ 
line  data  given  by  Fels  and  Schwazkopf  (1981).  In  Fig.  3,  infrared 
cooling  rates  due  to  the  15  ym  C02  band  are  shown  utilizing  standard 
and  polar  temperature  profiles  (Ou  and  Liou,  1983).  Also  included  in 
the  figure  are  a  number  of  published  results.  From  these  comparisons, 
it  is  noted  that  the  present  results  differ  from  those  of  Fels  and 
Schwarzkopf  (1981)  who  employed  a  line-by-line  technique  by  no  more 
than  5%.  It  should  be  noted  that  the  new  C02  emissivity  parameteriza¬ 
tion  will  affect  flux  and  cooling  rate  calculations  only  above  about 
35  km. 

Moreover,  we  present  total  cooling  rates  due  to  HgO,  C02  and  03 
computed  from  the  present  parameterization  program  and  compare  them 
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with  those  calculated  from  a  more  elaborate  band-by-band  method.  This 
comparison  is  shown  in  Fig.  2b.  Also  depicted  in  the  figure  is  the 
observed  cooling  profile  in  the  tropics  presented  by  Cox  (1969). 

Although  the  atmospheric  water  vapor  and  temperature  profiles  correspond¬ 
ing  to  the  observed  cooling  rate  may  not  be  exactly  the  same  as  those 
for  the  tropical  atmosphere  used  in  theoretical  calculations,  the  com¬ 
parison  does  serve  as  a  general  check  on  the  reliability  of  the  theoret¬ 
ical  result.  Results  from  the  broadband  parameterization  method  are 
less  than  the  band-by-band  values  by  ~0.2-0.5°C  day"1.  Moreover,  we 
see  that  the  theoretical  cooling  rates  are  in  general  agreement  with  the 
observed  values  in  the  troposphere.  In  the  stratosphere  the  total 
cooling  profile  derived  from  the  broadband  approach  shows  a  close  com¬ 
parison  with  that  from  the  band-by-band  method,  with  differences  gener¬ 
ally  less  than  about  0.2°C  day"1.  In  a  clear  tropical  atmosphere  the 
infrared  cooling  profile  shows  a  number  of  maxima  and  minima,  in 
addition  to  the  strong  cooling  in  the  upper  atmosphere  due  to  carbon 
dioxide  and  ozone.  The  minimum  at  ~17  km  is  known  as  a  basic  result 
of  the  temperature  inversion  at  the  tropopause.  The  maximum  at  -10  km 
and  the  minimum  at  -3  km  are  associated  with  the  nonlinear  interaction 
of  temperature  and  emissivity  slopes.  The  latter  slope  is  further 
associated  with  the  water  vapor  concentration  profile. 

In  addition,  in  Fig.  4  is  shown  an  infrared  cooling  rate  profile 
generated  from  the  present  broadband  emissivity  method  for  the  standard 
atmospheric  condition.  Also  illustrated  in  the  figure  for  comparison 
purposes  are  results  presented  by  Manabe  and  Strickler  (1964)  and  by 
Dickinson  (1973)  for  atmospheres  above  30  km.  Cooling  rates  derived 
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Comparisons  of  the  infrared  cooling  rates  and  net  fluxes  computed 
from  the  present  study  with  other  published  results.  Dots  in  the 
net  flux  diagram  denote  radiometersonde  observations  reported  in 
Elllngson  and  Gllle  (1978). 
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from  the  present  study  compare  well  with  Dickinson's  results.  We 
show  larger  cooling  rates  than  those  presented  by  Manabe  and  Strickler 
below  about  6  km.  Basically,  this  is  due  to  the  incorporation  of  water 
vapor  continuum  absorption  in  the  present  study.  The  general  differ¬ 
ences  between  the  present  study  and  that  of  Manabe  and  Strickler  may 
be  attributed  to  the  vertical  resolution  and  emissivity  method  employed 
in  the  calculation.  The  right-hand  graph  in  this  figure  illustrates 
the  infrared  net  flux  profiles  for  a  tropical  atmosphere  computed  from 
the  present  study.  Also  depicted  in  the  diagram  for  comparison  purposes 
are  results  presented  by  Ellingson  and  Gille  (1978).  The  dashed  line 
and  dots  are  net  flux  profiles  based  on  their  model  calculation  and 
radiometersonde  observations,  respectively,  reproduced  from  Fig.  14  of 

their  paper.  There  is  a  close  agreement  of  the  infared  net  fluxes 

„2 

between  model  computations  within  about  5  W  m  below  about  10  km. 

_2 

Between  about  10-24  km,  differences  increase  to  within  about  15  W  m  .  The 
present  results  computed  in  the  height  coordinate  appear  to  compare  more  closel 
with  observed  net  flux  data.  Note  that  the  computations  of  cooling 
rates  and  net  fluxes  utilize  two  different  atmospheric  profiles  so  as 
to  have  adequate  comparisons  with  other  published  results. 

Finally,  we  show  in  Fig.  5  infrared  cooling  rate  profiles  in  cirrus 
cloudy  atmospheres  employing  the  technique  described  in  Section  2.2 
and  compare  these  with  results  derived  from  a  band-by-band  model  (Liou 
and  Ou,  1981).  Tropical  atmospheres  containing  a  0.1  and  a  3  km  non¬ 
black  cirrus  with  a  base  height  of  8  km  were  used  in  the  calculations. 

For  the  0.1  km  cirrus  case  shown  in  Fig.  5a,  except  near  the  cloud  top 
and  bottom,  the  parameterization  approach  gives  reasonably  accurate 


cooling  values  as  compared  to  results  computed  from  a  more  detailed 
band-by-band  method.  The  larger  cooling  deviation  is  caused  by  the 
inability  of  the  broadband  calculation  to  obtain  very  accurate  net 
fluxes  near  the  thin  cloud  region.  The  most  significant  and  pronounced 
effect  of  non-black  cloud  is  the  suppression  of  cooling  below  the  cloud. 
The  cooling  rate  is  reduced  by  as  much  as  1°C  day"*  everywhere  below  a 
0.1  km  cloud.  The  cooling  profile  involving  a  3  km  cirrus  cloud  is 
shown  in  Fig.  5b  where  we  also  find  a  consistent  accuracy  of  the  broad¬ 
band  emissivity  approach  relative  to  the  band-by-band  method.  The 

_p 

vertical  ice  content  in  this  case  amounts  to  156.8  g  m  which  gener¬ 
ates  infrared  broadband  reflectivity  and  emissivity  of  0.2401  and  0.9876, 
respectively,  with  a  zero  transmissivity.  As  far  as  the  conventional 
definition  of  the  emissivity  is  concerned,  the  3  km  cirrus  almost 
approaches  a  blackbody.  However,  we  see  a  significant  component  of 
reflected  flux  from  an  almost  black  cloud  due  to  the  scattering  of 
infrared  radiation  by  cloud  particles.  The  contribution  of  reflection 
and  cloud  bottom  emission  produces  heating  instead  of  cooling  for  atmos¬ 
pheres  between  the  cloud  bottom  and  about  5  km  height.  The  reflection 
component  attributes  to  -20%  of  heating  in  this  case.  Moreover,  the 
surface  infrared  cooling  is  also  much  reduced,  by  ~2°C  day"*.  On  the 
other  hand,  atmospheric  cooling  is  enhanced  as  much  as  1-2°C  day"* 
between  20  and  30  km.  Above  30  km,  no  significant  variation  of  cooling 
rates  is  found  when  a  3  km  cirrus  is  present  in  the  atmosphere.  Further, 
for  internal  cooling  we  obtain  a  cooling  rate  of  -8°  C  day"*.  For  pur¬ 
poses  of  comparison,  we  also  depict  in  Fig.  5b  cooling  rates  derived 
from  the  band-by-band  method  employing  a  three-layer  cloud  model. 


Strong  cooling  (~10°C  day  )  at  the  cloud  top  and  rigorous  warming  (~6  C 
day’1)  near  the  cloud  base  are  especially  noticeable.  The  balance 
between  the  upward  and  downward  fluxes  near  the  cloud  center  results 
in  a  near  zero  cooling  rate.  The  broadband  cooling  rate  -8°C  day"1 
coupled  with  ~1°C  day"1  warming  at  the  cloud  base  appear  to  be  consis¬ 
tent  with  the  results  derived  from  the  band- by- band  method. 

Other  atmospheric  profiles  also  have  been  employed  in  performing 
broadband  and  band-by-band  infrared  cooling  calculations.  Similar 
agreements  between  the  parameterized  broadband  and  band- by- band  methods 
were  found.  Finally,  it  should  be  pointed  out  that  the  present  infrared 
parameterization  scheme  involving  nonblack  clouds  simplifies  the  intri¬ 
cate  cloud  interactions  with  atmospheric  infrared  fluxes  in  terms  of  the 
broadband  cloud  radiative  properties.  The  numerical  effort  and  computa¬ 
tional  time  requirement  are  consequently  equivalent  to  those  of  dear- 
column  cases. 
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Section  3 


SOLAR  RADIATION  PROGRAM 

3.1  Clear  Atmosphere 

In  a  clear  atmosphere,  we  consider  that  Rayleigh  scattering  due 
to  air  molecules  is  insignificant  above  a  height  Zj  of  10  km  for  the 
downward  solar  flux.  On  the  other  hand,  however,  the  upward  solar  flux 
above  this  height  is  contributed  from  the  Rayleigh  reflection  upward 
plus  a  diffuse  reflection  component  resulting  from  multiple  reflections 
between  the  Rayleigh  layer  and  the  underlying  surface.  Hence,  the 
downward  and  upward  solar  fluxes  at  level  z  above  Zj  may  be  expressed  by 

F+(z)  =  yQ  SQ  [1  -  A(zrz)]  ,  (3.1) 

Ff(z)  =  F+(Zl)  (ra  +  G  ta)  [1  -  A(z-Zl)]  ,  (3.2) 

where  z-j-  is  the  height  at  the  top  of  the  atmosphere,  the  solar  constant 
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Sq  is  taken  to  be  1353  W  m  ,  and  the  cosine  of  the  solar  zenith  angle 
is  defined  by 

u0  =  sine  sin6  +  cose  cos6  cosA  ,  (3.3) 

where  6  is  the  latitude,  6  denotes  the  inclination  angle,  which  Is  a 
function  of  the  day  of  the  year,  and  A  is  the  hour  angle.  The  flux 
reflection  r  and  transmission  t.  due  to  Rayleigh  scattering  is  taken 

a  a 

to  be  rg  =  0.28  /  (1  +  6.43  jjq)  (Lacis  and  Hansen,  1974)  and  t  *  1  - 
ra  and  ra  and  t.  are  the  corresponding  diffuse  (or  global)  values  which 

a  a  a 


are  to  be  obtained  from  integrating  y^  over  the  upper  hemi sphere.  The 
gaseous  broadband  solar  absorptivity  A(z)  will  be  defined  later  and  its 
diffuse  value  A(z)  =  A  [d  u(z)]  where  d  =  1.66  for  H^O  and  1.9  for  Og. 

The  higher  diffusivity  factor  d  for  ozone  is  chosen  to  account  for  Rayleigh 
scattering  effects  on  the  ozone  absorption.  Based  on  a  simple  adding  pro¬ 
cedure  the  multiple  reflections  involving  the  Rayleigh  scattering  layer 
and  the  surface  having  a  Labertian  albedo  of  r$  are  described  by  a  non- 
dimensional  upward  diffuse  component  given  by 

G  ■  ‘a  rs  (1  '  7a  rs  •  ■  <3'4> 

In  Eq.  (3.2),  the  ra  term  represents  the  direct  reflection  of  the  down- 

a 

ward  solar  flux  while  6  t_  denotes  the  diffuse  reflection  component. 

d 

Likewise,  the  downward  and  upward  fluxes  at  heights  below  Zj  may  be 
expressed  by 

F^z)  =  yQ  SQ  [1  -  A(z-j--z)]  ,  (3.5) 

F+(z)  =  F*(Zj)  [1  +  A(zj-z)]  ,  (3.6) 

It  should  be  noted  that  Rayleigh  scattering  below  Zj  will  reduce  both  up¬ 
ward  and  downward  fluxes  by  the  same  amount  since  its  overlap  with  gaseous 
absorption  in  the  near  infrared  may  be  neglected.  Thus,  the  net  flux  and 
hence  the  heating  rate  may  be  evaluated  without  the  incorporation  of 
Rayleigh  scattering. 

Three  absorbers  are  considered  in  the  transfer  of  solar  radiation, 
i.e.,  water  vapor,  carbon  dioxide  and  ozone  with  the  former  two  gases 
responsible  for  absorption  in  the  near  solar  infrared  region  and  the 
latter  in  the  UV  region.  It  is  assumed  that  overlap  of  HgO  and  C02 


absorption  and  0^  absorption  is  insignificant  so  that  the  broadband 
solar  absorptivity  may  be  expressed  by 

A(z)  =  |  A.(u 0/pQ)  f.  +  l  A.(Uw/y0)  f.  +  e  A5(uc/yQ)  ,  (3.7) 

where  uQ,  uw  and  uc  represent,  respectively,  path  lengths  for  03,  H20  and 
C02»  and  f..  are  the  fractional  solar  flux  associated  with  the  individual 
band  absorptivities  A^.  The  H20-C02  overlap  correction  is  given  by  e  * 
a  +  b  Ag(uw)  where  a  and  b  are  empirical  constants  derived  from  the 
measured  data  presented  by  Howard  et  al.  (1956);  a  =  0.75  and  b  =  -0.3/200 
for  A5(uw)  >  300  cm"1,  and  a  =  1  and  b  =  -0.7/300  for  Ag(uw)  <  300  cm-1. 

The  individual  absorptivity  for  ozone  is  given  by 

WM0^  =  1  "  exp  (‘ki  V'o*  *  (3.8) 

where  k^  denote  the  ozone  absorption  coefficients  taken  from  the  data 
obtained  by  Inn  and  Tanaka  (1953)  and  Vigroux  (1953)  anduQ  is  in  units  of 
cm  NTP"1.  These  coefficients  and  fractional  solar  fluxes  are  tabulated  in 
Table  3.  For  water  vapor  and  carbon  dioxide,  the  empirical  band  equation 
derived  by  Liou  and  Sasamori  (1975)  is  used.  It  has  the  following  analytic  form 


V^o'  ■  ^  [C1  +  Di  'o»10  (£  ^  +  X01]]  .  (3-9) 

where  Av..  are  the  band  widths,  and  Xgi  are  empirical  constants. 

These  constants  are  listed  in  Table  4.  The  pressure  P  in  units  of  mm  Hg 


is  weighted  over  the  path  length  given  by  P 


uc  is  in  units  of  g  cm 


-2 


fu 

=  P(u) 
Jn 


du/u  where  uw  or 


3 . 2  Cloudy  Atmosphere 

As  in  the  case  of  the  infrared  radiation  program,  we  consider  a 
general  two-layered  cloud  system.  We  wish  to  obtain  the  upward  and 
downward  solar  fluxes  at  levels  above  the  high  cloud,  between  the  high 


Table  3.  Ozone  Absorption  Coefficients  and 
Fractional  Solar  Fluxes 


AA(um) 

i 

O 

o 

i 

•r— 

fi  (Wm'2)* 

0.20-0.21 

1 

9.8 

1.24  x  10'4 

0.21-0.22 

2 

27 

2.97  x  10"4 

0.22-0.23 

3 

75 

4.59  x  10"4 

0.23-0.24 

4 

164 

4.59  x  10"4 

0.24-0.25 

5 

254 

5.14  x  10'4 

0.25-0.26 

6 

290 

7.55  x  10"4 

0.26-0.27 

7 

241 

1.35  x  10"3 

0.27-0.28 

8 

145 

1.59  x  10"3 

0.28-0.30 

9 

33.7 

6.46  x  10"3 

0.30-0.32 

10 

2.8 

1.01  x  10-2 

0.32-0.34 

11 

0.16 

1.50  x  10~2 

0.34-0.35 

12 

0.014 

7.95  x  10"3 

0.45-0.50 

13 

0.0107 

7.46  x  10"2 

0.50-0.55 

14 

0.055 

6.78  x  10-2 

0.55-0.60 

15 

0.11 

6.30  x  10-2 

0.60-0.65 

16 

0.09 

5.87  x  10"2 

0.65-0.70 

17 

0.038 

5.33  x  10‘2 

0.70-0.80 

18 

0.015 

9.14  x  10"2 

n 

*  Based  on  a  solar  constant  of  1353  W  m  and  the  fractional  fluxes 
presented  by  Thekaekara  (1974,  1976). 
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Table  4.  Empirical  Constants  for  1^0  and  CO^  Bands 
and  Fractional  Solar  Fluxes. 


A(ym) 

i 

C. 

D. 

1 

K. 

1 

Av..  (cm"*) 

1 

* 

CNJ 

1 

jj 

T* 

4- 

H^O  band 

0.94 

1 

-135 

230 

125 

1400 

0.0760 

1.1 

2 

-292 

345 

180 

1000 

0.0528 

1.38 

3 

202 

460 

198 

1500 

0.0732 

1.87 

4 

127 

232 

144 

1100 

0.0386 

2.7 

5 

337 

246 

150 

1000 

0.0242 

3.2 

6 

-144 

295 

151 

540 

0.0088 

C0o  band 

C 

2.7 

5 

-137 

77 

68 

320 

0.0073 

*  Based  on  a  solar  constant  of  1353  Wm"1"  and  the  fractional  fluxes 
presented  by  Thekaekara  (1974,  1976). 
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and  low  clouds  and  below  the  low  cloud.  Utilizing  the  approach  similar 
to  the  flux  transfer  in  a  clear  atmosphere  and  in  reference  to  Fig.  1, 
the  downward  and  upward  solar  fluxes  above  z^  may  be  written  in  the 


F^z)  =  pQ  S  [1  -  A (z)]  , 


Ff(z)  =  F4'(zu)  R3  [1  -  A(z  -  zfl)]  , 


(3.10) 


(3.11) 


z  >  z. 


where  R^  is  the  solar  flux  reflection  involving  the  two  cloud  system 

r, 

and  the  surface  defined  below  and  F  (z^)  is  the  downward  solar  flux 
at  position  ztl* 

To  derive  the  flux  from  below  the  high  cloud,  we  let  the  solar 
flux  reflection  and  transmission  for  the  high  (1)  and  low  (2)  clouds 
be  (r^,tj)  and  (r2,t2),  respectively,  with  their  corresponding  diffuse 
(or  global)  quantities  denoted  as  7j,  tj,  r2  and  t2  where 


f1  r<V 

' 0  t(p0) 


yo  %  * 


(3.12) 


Moreover,  we  also  let  and  A2  be  the  diffuse  gaseous  absorption  between 
the  high  and  low  clouds  and  between  the  low  cloud  and  the  surface.  We 
also  utilize  a  non-dimensional  total  upward  generation  function  G2^  which 
is  initiated  from  the  reflection  of  the  low  cloud  with  a  contribution 
from  the  surface  to  be  defined  below.  Thus,  the  downward  and  upward 
fluxes  between  zbl  and  zt2  may  be  written 

F*(z)  =  F^(ztl)  [tj  +  G^t}  (1  -  Aj)  7X]  [1  -  A(zbl  -  z)]  ,  (3.13) 


yryr 


In  Eq.  (3.13),  the  tj  term  represents  the  direct  transmission  of  the 
broadband  solar  flux  reaching  the  high  cloud  top,  while  the  following 
term  denotes  the  diffuse  component  which  is  generated  through  the 
generation  function  G^. 

Similarly,  if  we  let  the  direct  transmission  for  the  two  cloud 
system  be  r2  (the  corresponding  reflection  will  be  denoted  by  R2)  and 
the  upward  generation  function  originating  from  the  surface  be  G^, 
then  the  downward  and  upward  fluxes  below  the  low  cloud  base  are 
described  by 

F+(z)  =  F*(zt)  T3  [1  -  A(zb2  -  z)]  ,  (3.15) 

Ff(z)  =  F*(zt)  G3  [1  -  A(z)]  ,  (3.16) 

zb2  '  z  >  0 

where  T 3  denotes  the  direct  transmission  for  the  system  of  two  clouds 
and  the  surface. 

The  derivation  of  the  flux  reflection  and  transmission  and  the 
upward  generation  function  for  multiple  cloud  layers  is  similar  to  the 
adding  method  for  radiative  transfer  (see,  e.g.,  Liou,  1980;  Hansen, 
1971)  except  that  we  have  separated  the  direct  (associated  with  yQ) 
and  diffuse  components  for  the  transmitted  beam.  In  this  approach 
we  postulate  that  the  broadband  solar  flux  either  in  the  form  of  direct 
or  diffuse  beam  is  additive  by  means  of  multiple  reflections  between 
layers  concerned.  Thus,  the  flux  reflection  and  transmission  for  two 


cloud  layers  are 


r2  *  »*i  +  62  (1  -  Ax)  tj  ,  (3.17 

Tz  =  tj  (1  -  Ax)  t2  +  G2  (1  -  Aj)2  7j  t2  ,  (3.18 

where  the  first  and  second  terms  represent,  respectively,  the  direct 
and  diffuse  contribution  with  the  latter  resulted  from  the  upward 
generation  function  from  the  low  cloud  layer  given  by 

62  =  h  (1  "  r2  C1  -  0  -  V2  rx  r2]_1  ,  (3.19 

with  the  last  term  denoting  multiple  reflections  between  two  clouds. 

It  follows  that  the  reflection  and  transmission  for  a  combination  of 
two  clouds  and  the  surface  is  then  given  by 

R3  •  r2  +  63  (1  -  A2)  F*  ,  (3.20] 

r3  =  Tz  +  G3  (1  -  A2)  R*  ,  (3.21] 

where  the  upward  generation  function  from  the  surface  is 

G3  =  Tz  (1  -  A2)  r$  [1  -  (1  -  A2)2  R*  rs]_1  ,  (3.22] 

_ it  _ it 

and  R2  and  Tz  are  the  diffuse  reflection  and  transmission  for  the  two 
cloud  system  when  radiation  is  incident  from  below.  Thus,  Eqs.  (3.17) 
(3.19)  may  be  modified  to  yield 

R2  =  r 2  +  G2  (1  -  Ax)  t2  ,  (3.23] 

?*  -  \  (1  -  Aj)  tj  ♦  G*  (1  -  A,)2  r2  tj  , 
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(3.24 


(3.25) 


~  (1  -  A^)  [(1  -  A^)  Tj  Tg] 

Finally,  in  Eqs.  (3.13)  and  (3.14),  the  total  upward  generation  function 
contributed  from  the  low  cloud  and  the  surface  is 

G2^  =  62  +  G3  ^  *2  ^  "  ^1^  **1  ^2^-1  *  (3.26) 

where  the  second  term  is  required  in  order  to  give  proper  internal  fluxes. 

In  the  case  of  one  cloud  layer,  only  Eqs.  (3.10)  -  (3.14),  (3.17) 
and  (3.19)  are  applicable.  In  this  case,  the  total  diffuse  generation 
function  G^  in  Eqs.  (3.13)  and  (3.14)  should  be  replaced  by  the  value 
G2  defined  in  Eq.  (3.19)  and  the  total  reflection  at  the  top  of  the 
cloud  system  R3  should  be  replaced  by  the  value  R?  defined  in  Eq.  (3.17). 

The  broadband  solar  reflection  and  transmission  values  for  various 
cloud  types  are  obtained  through  the  following  parameterization  equation 
based  on  detailed  radiative  transfer  calculations: 

V>  ■  L  L  <?  W" 

m=u  n-u 

where  b  are  empirical  coefficients  derived  based  on  the  parameteriza- 
mn 

tion  principle  proposed  by  Liou  and  Wittman  (1979),  and  the  vertical 
liquid  content  W  is  in  units  of  102  g  m'2.  The  bmn  for  Cu,  As,  St  and 
Ci  are  respectively  given  in  Tables  3,  4,  5  and  6  of  that  paper.  Finally, 
computations  of  solar  heating  rates  in  clear,  overcast  and  partly  cloudy 
conditions  follow  Eqs.  (2.20)  and  (2.21). 

3 . 3  Verification  of  Solar  Parameterization  Program 

In  this  section,  we  present  some  results  calculated  from  the  solar 
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parameterization  program  described  in  previous  sections.  Comparisons 
will  also  be  made  with  results  computed  from  a  more  exact  and  compre¬ 
hensive  radiative  transfer  program  developed  at  the  University  of  Utah 
over  the  last  few  years  (Freeman  and  Liou,  1979;  Liou  and  Wittman,  1979). 

In  Fig.  6  (Liou  and  Ou,  1983)  are  shown  the  comparisons  of  the 
solar  net  fluxes  and  heating  rates  computed  from  the  parameterization 
scheme  and  a  more  exact  and  comprehensive  transfer  program.  A  fraction 
of  daytime  of  0.5  along  with  values  depicted  in  the  diagram  are  used 

in  the  computation.  Agreement  between  these  two  methods  is  within 

-1  -2 
~0.02°C  day  (solar  day  is  used  here)  and  5  W  m"  for  solar  heating 

rates  and  net  fluxes,  respectively.  Below  about  10  km,  solar  heating 

is  basically  produced  by  absorption  of  H20.  Between  about  10-15  km, 

absorption  of  C02  contributes  to  within  about  0.1°C  day"1  heating  rates. 

Above  15  km,  large  solar  heating  is  exclusively  caused  by  03  absorption. 

_2 

The  net  solar  flux  ranges  from  about  350  W  m  at  the  top  of  the  atmos- 

_2 

phere  to  about  270  W  m  at  the  surface.  A  series  of  comparisons  be¬ 
tween  the  parameterization  and  more  exact  programs  involving  different 
atmospheric  conditions  and  solar  zenith  angles  have  also  been  performed. 
We  found  that  the  present  parameterization  program  for  net  flux  and 
heating  rate  computations,  which  has  not  been  presented  in  our  previous 
publications,  has  an  accuracy  similar  to  the  case  depicted  here. 

Next,  we  compare  the  net  solar  flux  and  heating  rate  profiles  in 
cloudy  atmospheres  computed  from  the  present  parameterization  program 
with  those  obtained  from  a  more  exact  radiative  transfer  calculation. 

In  Fig.  7  (Hutchison  and  Liou,  1983)  a  1  km  single  cumulus  cloud  was 
inserted  into  a  midlatitude  winter  atmosphere  with  a  base  height  of 


2  km.  A  cosine  of  the  solar  zenith  angle  of  1,  a  surface  albedo  of  0.15 

_p 

and  a  solar  constant  of  1353  W  m  were  used  in  the  calculation.  The 
fractional  duration  of  sunlight  used  was  0.5  in  this  exercise.  Agree¬ 
ments  of  heating  rate  and  net  flux  profiles  between  the  parameterization 
method  and  a  more  exact  method  are  generally  excellent  in  the  cumulus 
case.  The  presence  of  a  low  cloud  clearly  alters  the  solar  heating 
profile  above  and  below  the  cloud.  A  heating  rate  of  more  than  10°C/24 
hr  is  produced  by  a  1  km  cumulus  cloud,  whereas  heating  rates  are  greatly 
suppressed  between  the  cloud  and  the  surface.  Similar  results  and  com¬ 
parisons  are  shown  for  the  middle  cloud  case  depicted  in  Fig.  8.  Because 
of  the  larger  liquid  water  content  in  the  altostratus  than  that  in  the 
cumulus,  a  larger  heating  rate  is  produced,  coupled  with  a  greater 
reduction  of  the  net  solar  flux  available  at  the  surface. 


altostratus 


Section  4 


PRELIMINARY  GLOBAL  RADIATION 
AND  CLOUD  RESULTS 


4.1  Global  Cloud  Cover 

In  order  to  test  the  applicability  of  radiation  programs  to  all 
atmospheric  conditions,  a  cloud  formation  scheme  was  used  In  connection 
with  the  AFGL  six-layer  general  circulation  model.  The  scheme  is  similar 
to  that  proposed  by  Liou  and  Zheng  (1984).  From  the  model  relative  humidity 
defined  as  h  =  q/qs(T),  where  q  and  T  are  specific  humidity  and  temperature, 
respectively,  criteria  for  cloud  formation  and  fractional  cloud  cover  are 
devised  on  the  basis  of  empirical -statistical  means.  When  the  predicted 
relative  humidity  h  is  greater  than  a  prescribed  value,  it  is  assumed 
that  the  water  and/or  ice  cloud  may  be  formed  with  the  cloud  cover  com¬ 
puted  from  an  empirical  equation,  whereas  if  h  is  smaller  than  or  equal 
to  the  prescribed  value,  no  cloud  is  allowed  to  form.  Thus, 


where 


h  >  h„ 
c 

h  <  hc 


(4.1) 


g(h)  =  {[h  -  hc(o)]  /  [1  -  hc(o)]}2  (4.2) 


and  the  critical  relative  humidity  as  a  function  of  o  is  given  by 


where  a  *  P/P*,  P*  the  surface  pressure,  *nd  the  an  are  empirical  coef¬ 
ficients  having  values  of  aQ  =  1,  aj  =  -0.218,  a2  =  -0.3196,  and  a3  * 
3.464.  In  reference  to  Fig.  9  and  subject  to  the  above  conditions, 
clouds  are  permitted  to  form  on  layers  a5,  and  a3-  In  cases  of 
multi-layered  clouds,  the  total  cloud  cover  is  calculated  by  means  of 
a  minimum  overlap  scheme,  i.e.,  the  largest  cloud  cover  obtained  will 
be  the  total  cloud  cover. 

Utilizing  the  initial  field  of  0000  Z,  15  January  1978,  consisting 
of  temperature  and  specific  humidity  provided  to  us  by  AFGL,  cloud  fields 
are  generated.  It  should  be  noted  that  no  prediction  is  performed  in 
this  exercise.  On  this  day,  IR  cloud  pictures  for  the  northern  hemi¬ 
sphere  from  the  N0AA  satellite  were  available.  Depicted  in  Fig.  10a 
is  an  IR  cloud  picture  covering  the  time  period  from  0756  Z,  14  January 
1978  to  0756  Z,  15  January  1978  for  vertifi cation  of  the  model  cloud 
generation.  Three  frontal  band  clouds  are  seen  in  the  Pacific  ocean. 

In  addition,  large  cloud  covers  are  also  observed  on  the  eastern  part 
of  the  United  States  and  Canada  and  over  the  Asian  continent.  Clouds 
associated  with  the  intertropical  convergence  zone  are  visible  in  the 
tropical  regions.  The  computed  cloud  covers  illustrated  in  Fig.  10b 
generally  agree  well  with  the  cloud  patterns  seen  from  the  satellite, 
except  in  the  tropical  regions  located  at  about  120°  W  and  60°  E  where 
overestimates  of  the  cloud  cover  are  seen.  In  this  presentation,  the 
minimum  cloud  cover  contour  line  is  20%  and  shaded  areas  represent 
cloud  covers  greater  than  50%.  This  good  comparison  between  the  com¬ 
puted  and  satellite  observed  IR  cloud  covers  is  quite  encouraging. 
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4.2  Global  Radiation  Budget 


Using  the  AFGLGCM  with  the  incorporation  of  UU  radiation  programs 

and  the  aforementioned  initial  data,  global  distributions  of  outgoing 

IR  and  reflected  solar  fluxes  at  the  model  top  were  generated  and  they 

are  illustrated  in  Figs.  11  and  12.  The  IR  and  solar  flux  contours  are 
_2 

in  units  of  W  m  .  In  Fig.  12,  regions  above  the  zero  line  denote  the 
solar  night.  For  comparison  purposes,  global  cloud  cover  distributions 
are  also  presented  (Fig.  13).  Generally,  lower  IR  fluxes  are  associated 
with  colder  cloud  fields.  On  the  other  hand,  cloud  fields  also  produce 
larger  reflected  solar  fluxes.  Since  these  results  are  preliminary  and 
since  the  prime  objective  of  these  illustrations  is  to  demonstrate  that 
the  present  radiati on/cl oud  programs  in  the  context  of  AFGLGCM  are 
functioning,  detailed  physical  interpretations  of  these  graphs  will  not 
be  given  here. 


OUTGOING  IR  FLUX 


Fig.  11.  The  computed  global  distributions  of  outgoing  IR 
AF6L  six-layer  spectral  model  on  15  January  1978. 


Fig.  12-  Computed  global  distributions  of  reflected  solar  flux  (in  units  of  W  m"  )  at  the  top  of 
the  AFGL  six-layer  spectral  model  on  15  January  1978. 


Computed  global  cloud  covers 


Section  5 


SUMMARY 

In  this  report  we  have  summarized  the  infrared  and  solar  radiation 
parameterization  programs  developed  at  the  University  of  Utah.  The 
infrared  radiation  parameterization  utilizes  the  broadband  emlssivity 
approach  in  which  the  emisslvlty  for  the  individual  absorption  band 
is  derived  either  from  band  models  (water  vapor  and  ozone)  or  line-by- 
line  data  (carbon  dioxide).  Me  have  also  presented  basic  flux  equa¬ 
tions  for  an  atmosphere  containing  high  and  low  clouds.  High  clouds 
are  treated  as  non-black  whose  broadband  emlssivity,  transmissivity 
and  reflectivity  are  parameterized  in  terms  of  the  vertical  ice  content. 
On  the  solar  side,  parameterization  has  been  carried  out  by  developing 
a  broadband  absorptivity  for  water  vapor,  carbon  dioxide  and  ozone 
absorption.  In  line  with  this,  cloud  reflection  and  transmission 
properties  are  parameterized  in  terms  of  the  liquid  water  content  and 
solar  zenith  angle.  Based  on  broadband  solar  radiative  properties  for 
gases  and  clouds,  we  have  devised  equations  for  the  computation  of 
upward  and  downward  fluxes  for  clear  and  one-  and  two-layered  cloudy 
atmospheres.  We  have  presented  in  some  detail  basic  parameterized 
equations  and  relevant  coefficients  and/or  constants  used  in  the 
current  version  of  the  radiation  code.  Accuracies  of  all  the  param¬ 
eterization  programs  have  been  verified  via  more  comprehensive  and 
exact  radiative  transfer  calculations  described  in  this  report. 


Moreover,  we  have  provided  a  comprehensive  documentation  of  the 
radiation  parameterization  programs  and  listings  of  all  the  relevant 
computer  formats  for  use  in  general  circulation  models.  Specifically, 
we  have  described  in  some  detail  the  link  and  interface  between  the 
UU  radiation-cloud  programs  and  the  AFGL  six-layer  spectral  model.  We 
have  also  presented  some  preliminary  results  of  the  global  distribu¬ 
tions  of  clouds  and  outgoing  IR  and  reflected  solar  fluxes  generated 
from  the  AFGL  six-layer  spectral  model  with  the  incorporation  of  UU 
radiation-cloud  parameterization  programs. 


Appendix 


•> 


h 


DOCUMENTATION  OF  RADIATION  PARAMETERIZATION 
PROGRAMS 

A.l  General  Description  of  the  Programs 

The  basic  principles  of  calculating  both  IR  and  solar  radiation 
fluxes  and  heating  rates  have  been  illustrated  in  Sections  2  and  3. 

In  this  section,  we  document  the  detailed  computer  routines  for  the 
IR  and  solar  parameterization  programs  for  use  in  general  circulation 
models. 

The  entire  radiation  program  is  composed  of  four  subroutines 
denoted  as  PRERAD,  IRRAD,  IRCLR  and  IRCLD.  The  subroutine  PRERAD 
serves  as  a  driver  program  which  connects  the  main  part  of  radiation 
calculations  with  other  GCM  programs.  It  interpolates  the  model - 
generated  temperature  and  humidity  profiles  from  the  main  part  of 
the  GCM  on  the  layers  designed  for  radiation  calculations.  The 
program  also  determines  the  cloud  position  along  with  the  cloud 
fraction  and  thickness.  The  computed  radiative  heating  rates  for 
radiation  layers  from  the  subroutine  IRRAD  are  then  weighted  by 
density  to  get  the  mean  radiative  heating  rate  for  each  GCM  layer. 

The  subroutine  IRRAD  is  the  main  program  for  radiation  calcula-  * 

tions.  It  is  composed  of  IR  and  solar  parts.  In  the  IR  part,  the 

v 

path  lengths  in  pressure  coordinates  are  first  computed  and  then  ad¬ 
justed  in  accordance  with  a  number  of  temperature  and  pressure  correction 
schemes  denoted  in  Section  2.  Next,  the  coefficients  of  the  third-order 
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polynomial  for  IR  emissivities  are  obtained  by  interpolation  between 
the  prescribed  reference  temperatures  listed  in  Table  1.  The  IR 
program  is  then  connected  to  a  large  computational  loop  in  which  the 
radiative  fluxes  for  each  cloud  configuration  are  calculated  based  on 
the  following  three  atmospheric  conditions: 

a.  For  a  clear  sky,  the  subroutine  IRCLR  is  called. 

b.  For  low  or  middle  clouds,  the  subroutine  IRCLR  is  called 
twice  for  regions  above  and  below  the  cloud. 

c.  For  a  high  (cirrus)  cloud,  the  subroutine  IRCLO  is  called 
once. 

Lastly,  in  the  IR  program,  the  total  heating  rate  profile  is  obtained 
from  the  net  flux  divergence  using  a  centered  finite  difference  scheme. 

In  the  solar  part,  the  path  lengths  of  each  model  layer  for  water 
vapor,  carbon  dioxide  and  ozone  are  first  computed.  These  layered 
path  lengths  are  then  substituted  into  equations  for  the  computations 
of  gaseous  absorptivities.  Again,  the  solar  program  is  also  connected 
to  a  large  computational  loop  in  which  the  solar  net  fluxes  for  each 
cloud  configuration  are  computed  according  to  the  following  atmospheric 
conditions: 

a.  For  a  clear  sky,  the  solar  fluxes  are  obtained  using  the 
previously  computed  absorptivities  with  Rayleigh  scattering 
corrections  given  in  Section  3. 

b.  For  cloudy  conditions,  cloud  reflectivity,  transmissivity 
and  absorptivity  are  first  computed  utilizing  formulas  pro¬ 
posed  by  Liou  and  Wittman  (1979).  The  net  solar  flux  for 
each  radiation  layer  is  computed  according  to  formulas  given 
in  Section  3. 


In  the  final  portion  of  the  solar  part,  the  solar  heating  rate  is  com¬ 
puted  from  the  net  flux  divergence  using  a  centered  finite  difference 
scheme. 


The  subroutine  IRCLR  computes  the  IR  fluxes  in  clear  regions. 
Thus,  it  can  be  used  to  compute  the  entire  clear  column  between  model 
top  and  surface,  the  clear  region  between  the  low  or  middle  cloud 
bottom  and  surface,  or  the  region  above  the  cloud.  The  subroutine 
IRCLD  computes  the  IR  fluxes  in  a  column  containing  only  one  layer  of 
a  non-black  cloud.  It  is  applicable  to  a  single  non-black  cloud  con¬ 
fined  either  between  the  surface  and  model  top  or  between  the  black 
cloud  top  and  model  top.  The  computational  procedure  is  similar  to 
that  described  in  IRCLR.  However,  the  flux  components  required  in 
IRCLD  are  four  in  both  the  above-cloud  and  below-cloud  regions.  These 
four  components  include  two  emitted  fluxes  in  upward  and  downward  dir¬ 
ections,  one  reflected  flux,  and  one  transmitted  flux.  The  advantage 
of  the  method  used  here  is  that  all  these  components  can  be  computed 
based  on  one  generalized  formula.  Consequently,  the  programming 
burden  is  greatly  reduced,  and  the  program  if  further  vectorized  can 
become  extremely  efficient. 


A.  2  Interfacing  with  AFGLGCM 

A  schematic  flow  diagram  illustrating  the  connection  between  the  AFGLGCM 
(Brenner  et  al . ,  1982)  and  radiation  programs  is  depicted  in  Fig.  14.  Two  sepa 
rate  subroutines  referred  to  as  ZENITH  and  CLOUD  are  also  connected  to  AFGLGCM. 
The  subroutine  ZENITH  computes  the  global  distribution  of  the  solar 
zenith  angle  while  CLOUD  computes  the  cloud  cover  distribution  in  the 
three  cloud-forming  model  layers.  Bascially,  the  major  modifications 
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are  all  done  in  the  subroutine  LALOOP  contained  in  AF6LGCM.  Other 
minor  modifications  include  (a)  storing  the  surface  albedo  data  in 
array  SALB  by  a  read  statement  in  the  main  program,  (b)  adding 
statements  in  subroutine  SPTOGP  to  compute  surface  pressures,  and 
(c)  inserting  print  statements  to  list  the  values  of  total  heating 
rate,  radiative  heating  rate,  temperature  humidity  and  surface  pres¬ 
sure  at  selected  grid  points  and  time  steps. 

In  the  common  block  CZALB,  SZ  stores  values  of  the  cosine  of  the 
solar  zenith  angle  and  SALB  stores  surface  albedo  values.  In  the 
common  block  STORE,  RADTEM  stores  global  radiative  heating  rate  data, 

ICL  stores  global  cloud  cover  information,  and  KFLUX  stores  data  for 
outgoing  IR  and  reflected  solar  fluxes  as  well  as  net  radiative  fluxes 
at  the  top  of  the  atmosphere.  In  the  common  block  RADPGR,  RIRAD  stores 
data  for  radiative  heating  rates  for  two  latitudinal  circles  described 
in  LALOOP.  Finally,  in  the  common  block  ZONAV,  YYY  stores  zonal ly 
averaged  total  and  radiative  heating  rates,  IYYY  stores  zonally  aver¬ 
aged  temperature  and  humidity  values,  and  IPPP  stores  zonally  averaged 
surface  pressure  values. 

The  basic  executing  procedures  for  the  new  radiation-AFGLGCM  program 
are  similar  to  those  for  the  original  AFGLGCM  program.  Also  note  that 
the  new  program  has  been  tested  and  verified  using  the  CRAY-1  computer 
in  NCAR.  Finally,  the  complete  listings  of  the  radiation  subroutines 
PRERAD,  IRRAD,  IRCLR  and  IRCLD  along  with  modified  subroutines  LALOOP, 
NLPROD,  and  SPTOGP  in  the  AFGLGCM  main  program  may  be  obtained  from 
K.N.  Liou  at  the  University  of  Utah. 
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